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ABSTRACT.  TK5~WOrfc  nf  thima-prereding  reports  on  the 
linearized  theory  of  vented  hydrofoils  is  compTel 
detailed  exposition  of  the  procedure  for  calculating  hydro¬ 
dynamic  coefficients.  The  basic  assumption  is  steady,  plane 
potential  flow  with  two  free  streamlines;  one  detaching  at  an 
arbitrary  chordwise  position  and  the  other  detaching  at  the 
trailing  edge.  The  single  cavity  terminates  at  an  arbitrary 
distance  downstream  of  the  trailing  edge.  The  calculation 
of  the  key  smooth- entry,  open- model  parameters  is  illus¬ 
trated  by  considering  a  symmetrical  and  a  cambered  pro¬ 
file,  each  vented  at  30$  chord.  The  theoretical  Cp  and  Cj_^ 
curves  for  two  vented  hydrofoils  with  exhaust  at  30$  chord 
are  compared  with  water- tunnel  test  data.  The  agreement 
is  not  as  good  as  in  previous  comparisons  for  leading-edge 
or  trailing- edge  free- streamline  detachment.  Real-fluid 
and  tunnel  blockage  effects  are  believed  to  have  produced 
the  discrepancies.  Because  the  calculations  for  arbitrary 
profile  involve  functions  that  are  independent  of  profile  form, 
the  function  tables  that  were  generated  during  this  work  are 
given  in  Part  Z  of  this  report,  bound  separately. 
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FOREWORD 


As  one  of  a  series  on  vented  hydrofoils,  this  report, 
consisting  of  two  parts,  completes  the  application  of 
linearized  theory  for  steady,  plane,  unbounded  flow.  The 
general  study  is  motivated  by  the  possible  use  of  gas  ex¬ 
haust  for  torpedo-control  purposes  and  the  use  of  base- 
vented  hydrofoils  as  propeller  blades. 

The  theoretical  calculations  are  illustrated  by  treating 
two  vented  hydrofoils  for  which  experimental  results  were 
given  in  Water-  Tunnel  Tests  of  Hydrofoils  With  Forced  Ven- 
tilation,  by  T.  CT  Lang,  D.  A.  Daybell,  and  K.  E.  Smith 
"(t'JAVOkD  Report  7008,  10  November  1959). 

The  comparison  of  theory  with  experiment  leads  to 
new  insight  into  the  nature  of  the  experimental  flows. 

Because  the  calculations  for  arbitrary  profile  involve 
functions  that  are  independent  of  profile  form,  the  func¬ 
tion  tables  that  were  generated  during  this  work  are  given 
in  Part  2  of  this  report,  bound  separately.  Thus  other  re¬ 
searchers  may  consider  further  hydrofoil  profiles  with 
minimal  labor.  However,  since  the  tables  are  of  restricted 
interest,  Part  2  is  being  sent  only  to  a  small  number  of 
those  receiving  Part  1.  Others  who  wish  to  receive  Part  2 
may  write  to  Commander,  U.  S.  Naval  Ordnance  Test  Sta¬ 
tion,  China  Lake,  Calif.,  Attn:  Code  7506. 

This  work  was  done  under  Bureau  of  Naval  Weapons 
Task  Assignments  RUAW-4E40 1  / 2 1  6  1 /R009-01 -003  and 
RUTO-  3 E- 000/ 2 1 6  1  /R009-01-003,  problem  assign¬ 
ment  401.  The  report  was  reviewed  for  technical  adequacy 
by  Dr.  Blaine  R.  Parkin  of  the  RAND  Corp. 


Released  by 
T.  E.  PHIPPS,  Head, 
Research  Department 
3  0  Ma  rch  19t)  3 


Under  authority  of 
Wm.  B.  McLEAN 
Technical  Director 
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NOTATION 

\lt  -  e  -  n/(I-  l)e 

A 

-  e  +  \/(l  -  l)e 

b 

(1  +  a)  CT 

ibi ,  &2  +  , , . 

Coefficients  in  expansion  of  w(z)  for  z — »<*> 

c  +  id 

loo  =  A*iX 

cn 

2  f1" 

c_  =  —  1  H(0)  cos  n0  d0,  n  =  1,  2,  . . . 

w  Jo 

C' 

Closure  coefficient 

cz 

z- plane  closure  singularity  coefficient 

CD>  cL,  cM 

Drag,  lift,  and  moment  coefficients  for  unit  chord 
length  and  unit  span  (Fig.  1 ) 

Cq 

Gas- supply  coefficients,  taken  from  Ref.  11 

dn 

cn  +  2(q  -  ag)  cos  n0' 

e 

Exhaust  location  measured  from  leading  edge  for 
unit  chord  length 

F,  G 

Functions  of  c,  d,  and  cos0  (page  12) 

1*  f2*  f3*  F4 

Functions  of  c,  d,  and  cos  0  (page  12) 

h 

Hydrofoil  contour  height  above  x  axis 

h',  h" 

Parameters  in  polynomial  approximation  of 
uniform- load  camber  line 

H 

Contour  slope,  dh±/dx 

i 

\TT 

In 

rQ'  1  cos  0  -  cos  0' \n 

-  d0 

JO  \  1  +  cos  0'  / 

iv 
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Im.  Imaginary  part  of 

T  it  /  cos  0  -  cos  O'  \  n 

Jn  -  dO 

J0'  \  1  +  cos  O'  / 

K  Cavity  number,  -  pc  >/pui 

f  Distance  from  leading  edge  to  cavity- termination  point, 
called  cavity  length 


^n  In  "  ^n 

Nu^  -  (a  Mb)  =  c  -  a  t  i(d  -  b) 

P  Vented  flat- plate  solution  for  unit  angle  of  attack  and  open- 
cavity  model 

Pn  In  t  Jn 

Q  Vented  flat- plate  solution  for  unit  angle  of  attack  and 
closed-cavity  model 

Rl.  Real  part  of 


tw  Wake  thickness  for  unit  chord 


I  tw  /  tW(j 

u,  v  x  and  y  components  of  velocity  perturbation  divided  by  Ujq 

w  Conjugate  complex  perturbation  velocity,  u  -  iv,  and 
analytic  function  of  t.  and  t, 

x,  y  Cartesian  coordinates  in  physical  plane  with  free  stream 
in  positive  x  direction  and  (0,  0)  at  leading  edge 

/.  x  -  iy 


n  Angle  of  attack 

Flat-plate  solution  parameter 

.  Parameters  in  polynomial  approximation  of  uniform- load 
camber  line 

f)  Vented  wedge  full  apex  angle 
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8-l/2*  80*  61*  •  •  • 


i 

e 

e' 


X..A 


v,  N 
P 


Parameters  in  expression  for  NACA  0010  profile 
slope  distribution 

Circle-plane  complex  coordinate.  £  =  &i  +  i t,i 
Circle- plane  polar  angle 

Value  of  6  at  leading- edge  point  on  unit  circle 

Coo  =  AeiX 

^-(a  +  ib)  =  Neiv 

Water  density 

Tunnel  height  divided  by  hydrofoil  chord 


W) 


Arctan  [  d  +  (-)Jl  -  a^  ]/(c 


a) 


Subscripts 

c  Cavity  value 

cl  Closed  model 

cw  Cavity- wall- speed  linearization 

fs  Free- stream- speed  linearization 
o  Open  model 

P  Open- model  flat- plate  solution  for  unit  angle  of 
attack 

Q  Closed- model  flat- plate  solution  for  unit  angle 
of  attack 

a  Smooth- entry  condition  for  open  model 

±  Respectively  upper  and  lower  sides  of  singularity 
slit  or  hydrofoil 

«  Condition  or  value  at  z  =  °° 


Superscripts 

+  Cusp  closure 

0  Zero  angle  of  attack 

—  Thickness  profile  and  its  parameters  after  re¬ 
moval  of  leading- edge  singularity 

~  Open- model  point- drag  profile  and  its  parameters 


vi 
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INTRODUCTION 


Theoretical  analyses  of  vented  hydrofoils  using  the  linearized 
"thin- airfoil"  theory  of  steady,  plane,  ideal  flow  are  presented  in 
Ref.  1  and  2.  This  third  report  presents  sample  theoretical  results, 
comparison  with  experiment,  and  function  tables  with  which  hydrody¬ 
namic  coefficients  can  be  calculated  according  to  the  formulas  devel¬ 
oped  in  Ref.  1  and  2. 

The  assumed  flow  configuration  is  shown  in  Fig.  1,  where  the  unit- 
chord  profile  has  a  full  cavity  (i.  e.,  f  >  1)  and  the  more  forward  free- 
streamline  detachment  point,  without  loss  of  generality,  is  assumed  to 
be  on  the  upper  side.  Thus  positive  lift  is  directed  away  from  the  foil 
on  the  vented  side,  and  positive  moment  tends  to  rotate  the  foil  so  as 
to  shield  the  vented  side  from  the  free  stream. 


FIG.  1.  Typical  Vented  Hydrofoil  Flow  Configuration  for  a 
Full  Cavity. 

The  tables  in  Part  2  of  this  report  are  in  four  groups.  Tables  1 
through  4  give  functions  convenient  for  treating  Cp,  C^,,  and  Cw  for 
a  broad  class  of  profiles  for  t  =  ».  Tables  5  through  8  give  the  func¬ 
tions  Fj,  F^,  F3,  and  F^  that  are  multiplied  by  a  profile's  wetted  con¬ 
tour  slope  distribution  and  then  integrated  to  obtain  the  profile's  key 
parameters  (neglecting  moment)  of  as,  uc  ,  tw  ,  and  Cls-  These  in¬ 
tegrations  are  for  nine  values  of  I  from  ^  to  1.1.  Table  9  gives  x(0) 
which  is  used  to  transfer  the  profile  slope  to  the  unit  circle  for  the  in¬ 
tegrations.  Auxiliary  tables  useful  in  the  calculations  and  in  possible 
future  work  are  Tables  10  through  14.  Tables  15  through  20  allow 
calculation  of  CD.  Cl.  and  K  versus  f  for  the  arbitrary  profile  at 
arbitrary  attack  angle. 
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In  this  work,  the  exhaust  point  location,  e,  is  in  effect  assumed  to 
be  in  the  range  0  <  e  <  1,  since  for  trailing-edge  venting,  e  =  1 ,  the 
work  in  Ref.  3  is  complete  in  itself.  However,  for  various  reasons, 
particularly  the  arrangement  of  tables,  the  case  of  e  s  1  is  also  in¬ 
cluded  here  at  several  points. 


INFINITE  CAVITY  LENGTH 


In  Ref.  2  the  relations  for  hydrodynamic  coefficients  for  I  *  »  are 
applied  mainly  to  profiles  of  simple  form.  The  function  tables  given 
here  allow  ready  treatment  of  a  fairly  broad  class  of  profiles. 

The  relations  for  1  =  oo  are 


CD  -  — 

2  ' 

(1) 

ir 

cM  -  — 

4 

-2  cos  0  '(3  +  4  cos2  0  ')di 


1  +  sJT\ 2  i 

-  di,  Cl  =  ir 


1  +  «/7\2 


(2di  cos  0'  -  d2) 


1  +  \Fe 


d4  -  6  cos  0'  d3  +  2(1  +  6  cos2  0')d2 


The  dn  coefficients  are  determined  by  the  contour  slope  and  the  angle 
of  attack  a,  with 

2  r  tr 

dn  =  cn  +  2(q  ”  os)  cos  n0\  cn  =  —  /  cos  n0  Hd0  (n  =  1,2,...) 

it  J  0 


(2) 


1  /**  1  -  Ve 

a8  =  _  I  H  d0,  cos  0'  - - -  a 

tt  ^  0  1  +  Ve" 


It  is  to  be  noted  that  as,  the  smooth- entry  attack  angle,  is  denoted  by 
a0  in  Ref.  2,  but  in  this  report  the  subscript  o  is  reserved  for  the  open- 
model  condition.  H(0)  is  the  wetted- surface  contour  slope  distribution, 
dh±/dx,  specified  over  the  upper  half  of  the  ^-plane  unit  circle.  The 
upper  surface  (+)  corresponds  to  0  <  0  <  0',  and  the  lower  surface  (-) 
corresponds  to  0'  <  0  <  ir,  according  to  the  wetted- surface  mapping  for 
f  =  *>, 


(3)  V7  =  ±| 


1  +  \J~e  1 


(cos  0  -  cos  0')  =  ± 


cos  0 

i  1  + 


-  cos  0'  \  /  cos  0  -  a  | 


cos  6  1 


\  1  +  a 


* 


See  Item  1, 


Appendix  A. 
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It  will  be  seen  later  that  these  uses  of  ag  and  a  for  I  =  oo  are  con¬ 
sistent  with  the  definitions  of  the  same  symbols  for  finite  I.  For  each 
profile  and  i,  a  -  a8  is  important  to  the  character  of  the  pressure  dis¬ 
tribution  near  the  leading  edge,  and  thus  ag  is  a  useful  parameter. 
Therefore,  for  continuity  between  the  results  for  infinite  and  finite  I, 
the  formulas  of  Eq.  1  are  put  in  the  form 


CD 


+  (a  -  ag)  V2rr 


1  -  \[e 
2 


z 

,  Cl  =  Cls  +  2tr(a  -  a8) 


1  +  \/T  2 
2  , 


(4) 


tr 

CM  =  Cms - (a  -  as) 


1  -  \Ze" 
1  +  >/e 


* 


The  smooth- entry- angle  parameters  are  seen  to  be 


(5) 


(6) 


(1T!)/o"cos9H 


de 


'1  +  -s/e\2  rir 


C,  =  2  — 


(7)  cMg  = 


1/1  +  N/e\4 


2  2 


/; ' 


a 


(Z  cos  0'  cos  0  -  cos  20)  Hd0 


cos  40-6  cos  9'  cos  30  +  2(1  +  6  cos^  0')cos  20 


-2  cos  0 ' (3  +  4  cos4  0')cos  0  ]  H  d0 


When  only  the  results  for  t  =  <x>  are  of  interest,  the  natural  formula¬ 
tion  in  terms  of  zero  angle  of  attack  is 


CD 


(8) 


=  yj4>+  Q^Tr|--^  ) 


1  +  \Te  \  z 


0  *  . 
CM  =  CM  -  -  Ql  “ 


1  +  \Te  \  4 


:L  =  CL  +  Zna  | 
/  1  - 


1+4 


1  +  \Te  , 


Noting  that  dn  =  cn  for  a  =  as,  and  using  the  relations 


* 


See  Item 


3, 


Appendix  A. 
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cos  29  -  cos  29'  =  4a  (cos  9  -  a)  +  2  (cos  9  -  a)2 
cos  39  -  cos  39'  =  (12a2  -  3)(cos  9  -  a)  +  12a  (cos  9  -  a)2  +  4{cos  9  -  a)3 
cos  49  -  cos  49'  =  (32a3  -  16a)(cos  9  -  a)  +  (48a2  -  8)(cos  9  -  a)2 

+  32a  (cos  9  -  a)^  +  8  (cos  9  -  a)^ 


obtain 

(9) 


Because  of  the  form  of  x(9)  in  Eq,  3,  any  term  in  H(9)  proportional  to 
(*/x)m,  where  m  is  an  integer,  leads  to  integrals  of  the  form 


Upper  side: 


Lower  side: 


/O  0 

The  integer  exponent  n  equals  m  +  1  inyCj},  m  f  1  or  m  f  2  in  Cj_,, 
and  m  +  1  through  m  f  4  in  Cjfa.  Tables  1  through  4  in  Part  2  give 
In.  Jn.  Pn  -  In  +  Jn«  and  Mn  ~  In  -  Jn  lor  m  -  l  through  10  and  for  9' 
values  corresponding  to  e  =  0,  0.  023,  0.  03,  0.1,  0.  2, ...  1.  0  plus  two 
additional  values  of  e  =  0.00516  and  0.  01754  corresponding  to  9’  =  30° 
and  40°.  The  evenly  spaced  values  are  the  same  as  those  used  in  the 
later  tables  for  finite  I.  The  two  extra  e  values  help  define  the  be¬ 
havior  for  e— *0,  which  can  be  very  nonlinear,  as  discussed  in  Ref.  2. 


See  Item  4,  Appendix  A 
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The  use  of  these  tables  will  be  illustrated  by  treating  typical  thick¬ 
ness  and  camber  distributions. 


NACA  00 10  THICKNESS  DISTRIBUTION 


The  first  example  is  the  NACA  four-digit  series  basic  thickness 
which  has  the  form  (Ref.  4) 

h±  =  ±5t(0.  2969  xi  -  0.  1260x  -  0.  35l6x2  +  0.  2843x2  -  0.  1015x4) 

where  t  is  the  maximum  thickness-to-chord  ratio.  Thus  for  the  NACA 
0010  profile, 

/  -1  2  3' 

—  =  ±(6.1  X-2  +  6q  +  6i  -  1  £-  •  1 


(12) 

where 


dx 


X  +  62  x^  +  63 


6.1  =  0.074225.  60  =  -0.06300,  =  -0.35160 

62  =  0.  42645,  63  =  0.  20300 

The  resultant  expressions  for  Cl- 

«  F  , 

(13) 


and  C^  are 


(14)  Cl  =  -4  +  6.  1  Pi  +  60M2  +  61M4  +  +  63M  gj 


(15) 


Cm  =  4 


(1  +  a)3 V  2 


8(1  +  a)2 


+  -  1  +  6qM^  +  6jMj  +  y  +  S^Mg) 

1  +  a  V  *  ‘ 

+  6_  *  P3  +  60M4  +  61Mb  +  6^Mg  +  63M10 

1 

Since  the  ±6- 1  x*  i  component  in  H  corresponds  to  the  point-drag  profile 
slope  for  1  =  uo,  there  should  be  no  contribution  from  this  term  to  C^, 
and  Cjjj.  This  is  seen  to  be  true  using,  in  Eq.  14, 
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a  air  a 

Pi  +  — —  Po  =  -  -  + -  (if)  =  0 

1+a  1+a  1+a 


and  in  Eq.  15, 


(1  +  a)- 


PO  +  -  p2  +  p3  s 


(1  +  a)J  1  +  a  l  (1+a) 


a  [ir(l/2- a^)]  air(3/2  +  a^) 


(1  +  a)“ 


Thus,  simplified  forms  for  the  0010  profile  are 

r~o  fi 

(13')  ^  Cp  =  /  —  ^ir  6.  i  +  6qMj  +  S^Mj  + 

(14')  C°  =  -4  -1-  |  ^0  _  ^  \  +  6qM2  +  6lM4  +  6zM6  +  63Mg 


Cm  =  4 


1  0 

ir6.A  + - Cl 

/  8(1  +  a)2 


+ - (^0^3  +  61M5  +  62M7  +  63M9) 

+  60M4  +  61M6  +  62M8  +  53M10 


Figure  2  gives  the  resultant  curves  of  Jc p,  C^.  and  Cj^  versus  e,  ob¬ 
tained  via  Tables  3  and  4.  For  another  thickness  ratio,  t\ 


;  resultant  curves  of  %/C]>  ^L’ 
3  and  4.  For  another  thickness 

-  10  *'  (\Z^)0010 


and  the  analogous  forms  for  and  cj[^. 

POLYNOMIAL  AND  UNIFORM-  LOAD  CAMBER  LINES 

The  second  example  is  the  fourth- order  polynomial  camber  line, 


Vl  ,  V2  ,  Y3  4  dh±  ,  3 

h±  =  Yo  x  +  —  x  +  —  x  +  —  x  >  -  =  Y0  +  Yl  x  +  Y2  x  +  Y3  x 

2  3  4  dx 


and  thus  (not  considering  moment), 
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In  Ref.  2,  this  type  of  camber  line  was  used  to  approximate  the  uniform 
load  or  NACA  "a  =  1.  0"  mean  line  (Ref.  4)  by  fitting  the  polynomial  to 
the  mean  line  for  design  of  1.0  to  give 


=  0, 

O 

If 

-H 

x: 

=■  0. 

25, 

=  0.  04475  =  h' 

=  0. 

50, 

=  0.  05515  =  h' 

=  0. 

75, 

=  0.  04475  ~  h' 

=  1. 

00, 

=  0 

The  five-point  fitting  produces 
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64 

Yo  =  -12  h'  +  — h".  Y1 
3 


640 

152  h' - h" 

3 


2 

V2  =  -384  h*  +  512  h".  Y3=  ’  -  Y2 


Figures  3  and  4  give  yCp  and  versus  e  for  the  fourth- order 
polynomial  approximation  to  the  uniform- load  camber  line.  Also 
shown  are  more  accurate  computer  solutions  for  the  precise  uniform¬ 
load  camber  line,  again  for  design  Cl  of  1.0.  These  were  obtained  by 

machine  calculation,  using  the  y  Cp  and  formulas  expressed  in 
terms  of  x  instead  of  8: 


*  \Tlir  (Jo  dx 


dx 


(18) 


r  1 .  . -  \  dh_ 

-  I  *fe  +  (1  -  \fe)  */x  -  x  I  - 

JO  1  ‘  dx 


dx 


r°  - 
CL  ~ 


(19) 


-  2  | - 1 +  J  's/x^'n/T  -  (1  -  *fe)  N fx  -  x  J 


-Idh^. 

dx 


dx 


f1  r  »  dh_  | 

I  \fx  l\fe  +  (1  -  •Je)  \Tx  -  xl  2 - dx  | 

Jo  1  dx  j 


The  machine  calculations  were  made*  on  a  Royal  McBee  IXrP-30  digital 
computer  by  using  uniform  x  increments  of  0.  005  and  taking  the  inte¬ 
grand's  value  at  the  midpoint  of  each  increment  as  its  average.  An 
important  increase  in  accuracy  was  obtained  by  analytically  correcting 
the  increments  adjoining  the  four  singular  points  at  x  =  0  and  e  in  the 
upper- surface  integrals  and  at  x  =  0  and  1  in  the  lower- surface  inte¬ 
grals.  To  calculate  in  the  same  way, 


1 


By  J.  A.  Poore  of  this  Station. 
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0.0  0.1  0.2  0.9  0  4  0.9  0.0  0.2  0.0  0.9  1.0 


FIG.  3.  Infinite -Cavity  Drag  Functions  for  Uniform- 
L,oad  Camber  Line  and  for  Approximating  Fourth- 
Order  Polynomial  Camber  Line. 


FIG.  4.  Infinite- Cavity  Lift  Functions  for  Uniform- 
Load  Camber  Line  and  for  Approximating  Fourth- 
Order  Polynomial  Camber  Line. 
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rt 


|<i  -  + 

+  J  x(l  -  sfe  +  2*/x)  [*/e  -  (1  -  •Je)\fx  -  x] 

-  J  x(l  -  <JT -  2 \GT)[*/e  +  (1  -  */e)  */x  -  x] 


id  K 


HYDROFOIL  B 

The  computer  calculation  was  also  used2  to  consider  the  NACA 
65^-015  profile,  with  modified  trailing  edge,  called  Hydrofoil  B  in 
Ref.  5.  Figures  5  and  6  show  the  resultant  \C £>  and  Cj^  for  Hydro¬ 
foil  B  and  for  the  approximation  made  to  it  in  Ref.  2  by  means  of  a 
roughly  similar  profile  with  a  very  simple  contour  expression. 


THE  SMOOTH-ENTRY,  OPEN- MODEL  PARAMETERS 

As  derived  and  illustrated  in  Ref.  1,  the  treatment  of  an  arbitrary 
profile  for  1  <  1  <  <®  reduces  in  the  main  to  calculating,  for  a  spread  of 
i  values,  the  several  parameters  of  the  smooth-entry,  open-model 
solution.  When  moment  coefficient  is  not  considered,  four  parameters 
are  involved,  whose  formulas  are 


*w8  .  ^Ls 

2x  4it  N 


•tu,=:lo’[AC 

■  1  r  [ 

N  |  w  JO  I 


AG  +  i 


F  Hd0 


/A2-l\  l2 
AG  +  i  I - If  Hde 


[2A^  sin  2\  +  i  (A*  -  1)] 
(A4  +  1  -  2A2  cos  Z\) 


(uca  +  i°s) 


where 


2  By  D.  A.  Daybell  of  this  Station. 
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FIG.  6.  Infinite- Cavity  Lift  Functions  for  Hydro 
foil  B  and  for  Similar  Profile  of  Simple  Contour 
Expres  sion. 
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1  1 

F  - - + - 

A2  +  1  .  i  A  cos  (\  +  0)  A2  +  1  -  2A  cot  (X.  -  0) 


2(A2  +  1)  -  4 A  cos  X  cos  0 


(A2-  1)^  -  4  ( A^+  1)A  cos  X  cos  0  +  4A*(cos^X  +  cos^  0) 

and 

sin  (X  +  0)  sin  (X  -  0) 


A^  +  i  .  2  A  cos  (X  +  0)  A2  +  1  -  2  A  cos  (X  -  0) 

sinX  [2  ( A2+  1)  cos  0  -  4Acos  x] 

(A2  -  l)2  -  4  ( A2  +  1)  Acos  X  cos  0  +  4A2  (cos^X  +  cos20) 


A  desirable  feature  of  these  integrals  is  that  there  are  no  singulari¬ 
ties  in  the  integrand  factors  of  H  as  long  as  1  <  f  £  00.  Thus  if  there  are 
no  singularities  in  H,  or  if  they  are  removed  as  described  in  Ref.  1,  the 
integrals  can  be  evaluated  by  simple  methods.  The  computations  for 
this  report  were  performed*  with  an  Ott  rolling-disk  planimeter.  For 
more  extensive  studies,  the  calculation  procedure  given  here  could 
easily  be  adapted  to  machine  computations. 

During  the  course  of  the  planimeter  work,  it  was  found  worth  while 
to  shift  to  the  following  format: 


(22) 


(23) 


"uc8  + 


r 

iwo8  = 

JO 


(F!  +  iF,?)  Hd0 


^Lg  hi  /*  TT 

tw„  + -  =—  e’lV  (F3  +  iF4)  HdG 

2  N  Jo 


with 


Fj  +  iF^  -  AG  +  i 


F3  +  iF^  -  (F i  +  iF^) 


■(¥) 


F1  +  iF2  •  .4 


2  sin  2  X  +  i  (A*  -  1) 


A’  +  1  -  2A4  cos  2X 


By  J.  M.  Hamilton  of  this  Station. 
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=  Fj*  .  F  £  + 


(A4-  1)F2  -  2 A2  sin  2^ 


+  i 


A4  +  1  -  2A2  cos  Z\ 

2A2sin  Z\  Fz  -  (A4-  1)  F* 


2Fi  F^  + 


A4  +  1  -  2 A2  cos  Z\ 


The  use  of  F3  +  iF4  avoids  some  unnecessary  compounding  of  errors, 
caused  by  subtracting  nearly  equal,  relatively  large  terms,  particularly 
for  large  <. 

The  planimeter  work  has  considered  the  following  spread  of  I 
values: 

00,  8,  4,  2.5,  2,  1.5,  1.25,  1.125,  and  1.1 

It  will  be  seen  that  the  absence  of  1.0  from  this  list  is  regrettable, 
since  the  s-parameter  curves,  when  plotted  versus  reciprocal  cavity 
length,  often  show  increasing  curvature  for  f — *1,  The  mapping  used 
caused  the  forms  of  the  s-parameter  integrals  to  be  indeterminate  at 
f  -  1,  but  it  can  be  seen  that  the  s-parameters  are  finite  there  unless 
the  profile  has  a  trailing- edge  singularity.  A  case  in  which  the  solu¬ 
tion  for  I  =  1  has  been  obtained  is  the  vented  wedge  discussed  in  Ap¬ 
pendix  B.  For  more  extensive  work,  it  would  be  useful  to  specifically 
treat!  ~  1  and  to  routinely  obtain  the  smooth-entry,  open-model  rela¬ 
tions  for  that  limit.  The  1  value  of  1.125  could  be  dropped,  keeping  the 
number  of  cases  to  nine. 


At  the  other  limit,  t  =  simplified  relations  are  already  available 
from  the  previous  work.  However,  it  is  necessary  to  confirm  that 
those  previous  results  are  recovered  with  the  more  general  forms. 
Using  the  notation 

Ae**'  -  c  +  id  -  (a  +  N  cos  v)  +  i  (b  +  N  sin  v) 


and  the  definitions 

v'f  -  e  -  \!  (i  -  1  )e 

a  = - - -  ,  b  =  (1  +  aWf  -  1 

\t  -  e  +  \/  (I  -  l)e 

which  lead  to 

_____  _____  I 

t  -  1  n/7  -  e  +  \/  e 

f 

it  can  be  seen  that,  for  i — ►*), 


N  a  (1  +  a)  4T\ 


l  -  1 


!  -  e 


13 


t«faW  *«t>HKi t*W»V.«*.  *  s 


NAVWEPS  REPORT  7941 
Part  1 


(24) 


1  -  </e  ’ 


c  *  2 


+  0(l_1),  d 


— -~W[l  +  0(l-1)] 
\  1  +  \Te  I 


i  +  *>/T  i 

The  expressions  for  Fj,  . . .  F4  in  terms  of  c,  d,  and  cos  6  become 

d[2(c2  +  d2  +  1)  cosG  -  4c] 


Fi  = 


F2  = 


<c2  +  d2  -  l)2  -  4{c2  +  d2  +  l)c  cos0  +  4c2  +  4(c2  +  d2)  cos20 

(c2  +  d^  -  l)(c^  +  d^  +  1  -  2c  cos0) 

(c2  +  d2  -l)2  -  4(c2  +  d2  +  l)c  cos0  +  4c2  +  4(c2  +  d2)  cos20 


,  [<c2  +  d2)2  -  1  ]  F2  -  4dcF i 

?3  =  Fi  -  F22  +  -- 


(c2  +  d2  +  l)2 


4c 


F4=2FlF2-  - 


[<c2  +  d2)2  -  l]  Fj  +  4dcF2 


(c2  +  d2  +  l)2  -  4c2 


These  can  be  expanded  in  powers  of  d,  with  the  coefficients  being  0(1) 
functions  of  I  and  0  for  I—**.  Thus  it  is  found  that 

Ft  =  (2  cos0)d_l  +  0(d-3) 

F2  =  1  +  2(1  +  c  cos0  -  2  cos2  0) d*2  +  0(d-4) 

F3  =  2(4  cos20  -  2  -  c  cos  0)  d-2  +  0(d'4) 

F4  =  (2  cos  0)d*1  +  0(d'3) 

Using  these,  and  also 

/  1  -  Ve 


cos  v 


j  — [l  +  Ofl’1)],  sinvrl  +  Otr1) 
/  Vf 


then  Eq.  22  and  23  yield 


1  + 


(25) 


irucB  +  i^as  = 


Ve  \  1  r-n 

- cos0  [l  +  0(1- M]H  d0 

+  i  J  [l  +  0(1"  ^JH  d0 
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.  C^s 
(Z6)  tw  +  i - 


r 

=  i[l+0(£-1)]  [ (cos  ^  F 3  +  sin  v  F4)  +  i(cos  v  F4  -  sinvF3)]HdG 

Jo 

/*  irf  / 1  +  \fe  \ 

=  Lt  - 

Jo  l  \  2  I 


1  +  \ fe  \2, 


r  w 


+  i 


cos  0  +  0(1”  *) 
1  -  *fe  \ 


1  +  sTei 


c  cos  0 

cos  0-2  cos^0  +  1  + - +  0(1”  ) 


Hd0 


Thus  for  i  =  00  Eq.  25  is  seen  to  give  Uc8  =  0  and  as  in  accord  with  Eq.  2. 
Also  using  the  expansion  of  c  from  Eq.  24, 


in  accord  with  Eq.  6.  Similarly  one  obtains 


1  =  00 


2 

lim  (2u  ^  )  =- 

- »  JO  s  s  ^ 


cos  0  Hd0 


in  accord  with  Eq.  5.  It  can  be  seen  that  the  sign  of  V^Dg.  important  to  the 
Cd  versus  q  variation  in  Eq.  4,  should  be  taken  the  same  as  that  of  tw  in 
Eq.  26,  so  that 


(27) 


cos  0  HdO 


To  illustrate  the  characters  of  F],  ...  Fq,  Fig.  7  gives  their  curves  for 
e  =  0.  3  and  for  t  =  4,  2,  1.  25,  and  1.1.  Since  the  F's  are  fairly  simple 

functions  of  0,  their  values  were  first  calculated  at  10-degree  intervals  in 
0  <  0  <  180°  and  then  at  intermediate  values  near  critical  points.  Additional 
0  values  were  also  needed,  particularly  for  0  near  0\  to  sufficiently  define 
the  products  FjH,  ...  F4H.  Tables  5  through  9  give  the  10-degree  values  of 
the  F's  and  of  x(0)  for  e  -  0.  3,  0,  4,  and  0.  5  and  for  i  =  *>,  8,  4,  2.  5,  2,  1.5, 

1.  25,  1.  125,  and  1.1. 


To  aid  in  the  calculation  of  other  tables  of  the  F's,  Tables  10  and  11  give 
the  major  parameters  c  and  d  for  all  the  e  and  i  values  of  the  flat-plate 
parameter  tables  given  later.  Tables  12  and  13  give  parameters  a  and  b 
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which  enter  into  the  wetted  surface  mapping  relation  for  x(6)  for  ar¬ 
bitrary  I,  _ 

x  1 

(28)  -  =  ±—  (cos  9  -  a) 

-  x  b 


Table  13  gives  hl/N2  that  is  used  in  obtaining  tWg  +  iCj^/2  since 

bi  bi 

—  e_iv  =—  [<c  -  a)  -  i(d  -  b)l 
N  N2 


Typical  steps  in  the  treatment  of  the  NACA  0010  thickness  distribu 
tion  and  the  uniform- load  camber  line  are  described  below. 


NACA  0010  THICKNESS  DISTRIBUTION 


The  first  step  in  treating  a  profile  with  a  round  leading  edge  is  the 
removal  of  the  leading-edge  singularity  in  H(6).  As  discussed  in  Ref.  1, 
this  is  done  by  subtracting  the  appropriate  point-drag  profile.  From 
£q.  12  of  this  report,  it  pan  be  seen  that  the  singular  term  in  the  0010 
profile's  H(x)  is  ±6_i  x"2  with  6_J  =  0.  074225.  From  the  open-model 
point-draa  profile  expression  in  Ref.  1,  the  corresponding  slope  dis¬ 
tribution  H(x)  for  arbitrary  f  is 


U9) 


H(x) 


■*s FJ 


X 

1  -  - 

f 


Thus  for  the  0010  profile,  uc  =  0.  074225/n/7~ and  H  in  the  s-parameter 
integrals  is  replaced  by 


TT=  H  -  H  =  ± 


+  6q  +  6j  x  +  6^  x 


2 


+  63  x3 


Figure  8  show JTT|  versus  x  for  several  I  values.  Note  that  H  =  ±6q  at 
x  =  0,  so  that  H(0)  will  be  discontinuous  at  0  =  0'.  As  an  example  of 
the  corresponding  TT(0),  Fig.  9  gives  H(0)  for  e  =  0.  3  and  i  =  4.  For 
the  same  e  and  I,  Fig.  10  gives  the  corresponding  curves  of  FjTT,  .  .  . 
F4TT and  the  planimeter-determined  areas. 

When  the  resultant  s-parameter  values  are  identified  by  bars  and 
the  point-drag  contributions  by  tildes,  the  final  s- parameters  are  ob¬ 
tained  using 


(30) 


uCs  =  uCs  +  UcB,  Qg  -  as  +  ‘3’s  =  <*s 


tw8  =  tWj|  +TWg  =  TWs  +  trucf, 
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1.  6 


(a)  FX 


(c)  F3 


FIG.  7.  Integrand  Factors  in 
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(d)  F4 

s- Parameter  Integrals,  e  =  0.  3  for  Various  I  Values. 
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FIG.  8.  H(x)  for  NACA  0010  Profile  for 
Various  i  Values. 


a  jo  to  »*  no  i« 

e,  d«g 


FIG.  9.  H(0)  for  NACA  0010  Profile, 
e  =  0.3  and  i  =  4. 


FIG.  10.  8-Parameter  Integrands  for  the  H(0)  of  the  NACA  0010  Profile 


0.  3  and  f  =  4, 
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Figures  il  and  12  give  the  curves  of  the  s-parameters  corresponding 
to  H  and  H,  showing  how  important  are  the  tilde  contributions  to  uc> 
and  ty,-.  "Cpp  and  Cp.  have  been  evaluated  for  f  =  *>  using  Eq.  27 
to  check  continuity  with  the  values  for  finite  i. 

One  check  on  the  accuracy  of  the  calculations  can  be  made  as  fol¬ 
lows.  For  e  =  0.  3  and  1  =  the  analytically  determined  values  for  the 
0010  profile,  from  Fig.  2,  are 


cj)  =  -0.0059,  =  -0.306 


The  planimeter  results  are 


as  =  0.0414,  Cd8  =  (0.  0161)  ,  CL  =-0.163 


which  imply 


V4  =  Zit  a 8  | 


1  -  \Te  ' 


=  -0. 0068 


Cl  =  CL  -  2 it  Qg 


1  +  \Fe 


=  -0.  319 


/  0  0 

The  effect  of  the  errors  in  \Cj}  and  Cl  can  be  seen  to  be  small, 
as  follows.  Cj)  and  Cl  from  Eq.  8  can  be  put  in  the  forms 


Cjj  =  2ir  (a  -  a+)^ 


1  -  \/7\2 


a+  = 


1  +  \l~e  I 


Cl  =  -o  n)  — 


■Cl 


2tt 


1  +\fey 


where  ajo  is  a  for  Cl  =  0.  Thus  a+  is  in  error  by  about  +0.09  degree 
and  q|q  by  about  +0.18  degree.  Thus  the  error  due  to  the  planimeter 
calculation  inaccuracy  will  be  small  away  from  the  zero  drag  and  zero 
lift  angles. 
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UNIFORM- LOAD  CAMBER  LINE 

The  uniform- load  camber  line  has  slope  singularities  at  the  leading 
and  trailing  edges.  For  the  planimeter  calculations,  the  singularities 
were  eliminated  in  a  practical  and  realistic  way  by  making  H(x)  uniform 
in  the  intervals  0  <  x  <  0.  005  and  0.  995  <  x  <.  1,  and  by  choosing  the  H 
values  therein  to  match  the  exact  profile  contour  height  at  the  interval 
end  points.  Since  H(x)  is  antisymmetrical  about  x  =  0.  5,  Fig.  13  gives 
the  approximating  H{x)  only  for  0  <  x  <.  0.  5.  Note  that  for  e  =  0.  3  the 
discontinuity  at  x  =  0.  005  in  Fig.  13  will  appear  at  three  6  values  in  the 
s-parameter  integrands,  corresponding  to  x  =  0.  005,  y  =  0±,  and 
x  =  0.  995,  y  a  0- .  Figure  14  gives  the  resultant  curve  of  H(0)  for 
e  =  0.  3,  t  -  4.  There  is  no  sign  reversal  at  9  =  0'  for  a  camber  line. 
Figure  15  gives  the  resultant  curves  of  FjH,  . . .  F^H  and  the  corre¬ 
sponding  areas,  which  lead  to  the  s-parameter  values  for  1/i  =  1/4 
seen  in  the  final  curves  in  Fig.  16. 

In  order  to  check  the  combined  accuracy  of  the  H  approximation 
and  the  planimeter  integration,  the  results  for  I  =  *  are  again  com¬ 
pared.  The  accurate  machine  computations  for  e  =  0.  3  shown  in  Fig.  3 
and  4  give 

=  0.  142,  C°  =  0.  324 
The  planimeter  results  for  f  =  <»are 

a*  =  0.0799,  ^Cd8  =  0.  1815,  CLs  =  0.  634 
and  these  lead  to 

=  0.  136,  C°L  =  0.  333 

Proceeding  as  for  the  0010  profile,  it  is  found  that  the  a+ error  is  about 
+0.  66  degree  and  that  of  a  for  =  0  is  about  -0.14  degree. 

While  the  accuracy  of  the  planimeter  integrations  has  been  shown 
in  two  cases  for  i  =  «o,  it  is  necessary  also  to  check  on  the  integrals 
and  the  planimeter  calculation  accuracy  for  finite  i.  This  has  been 
done  in  Appendix  B  by  considering  the  vented  wedge  and  comparing  its 
s-parameter  values  from  planimeter  integration  with  those  from  the 
formal  solution  constructed  of  point  singularities  on  the  unit  circle. 
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For  each  case  but  the  first, 


=  2ucp  twp ,  and  for  I  =  », 
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FIG.  15.  8- Parameter  Integrands  for 


of  Uniform- Load  Camber  Line,  e  =  0.3  and  1-4 
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FIG.  16.  Planimeter  Results  for  s- Parameters  Ver¬ 
sus  1/1  for  Uniform- Load  Camber  Line,  e  =  0.  3 

CDp  =  |2ucpVp  J  =  2"  |  “  J 

In  the  last  case,  e  =  1,  the  solutions  for  all  i  have  collapsed  to  the 
solution  for  fully  wetted  flow. 


CLOSED  MODEL 

The  key  parameters  of  the  closed-model,  flat-plate  solution  Q(C), 
also  for  unit  angle  of  attack,  are  given  in  Ref.  1  as 


u 


CQ 


Pq 


ucp 

Pp 


+  C'  (  A  +  — |  sin  \ 


1 


1 


•  - =  —  -  C'  I  A  -  — 

Pq  Pp  \  A 


(32)  CDq  =  2ttC2  =  8ir(bC'PQ)2f, 


cLq=~  cLp-  PQC' 

Pp 


1 


sin(\  -  v)  -  — —  sin(\  +  v) 
A 2 


blA 

2N 


where  the  closure  condition  gives 
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(53) 


C' 


1 

cos  (\  -  v)  +  — -  cos  (V  +  v) 
A4 


Tables  18,  19,  and  20  give  uCQ,  Cdq,  and  for  the  same  e  and  t 
values  as  the  P-parameters. 

For  two  of  the  four  boundary  cases,  the  closed  and  open  models 
are  identical,  because  the  closure  singularity  is  unimportant  for  i  =  °° 
and  is  zero  for  e  =  1.  For  i— *1,  the  closure  singularity  increases 
without  limit  and  Ucq  and  Clq — »*>.  The  Cdq  behavior  for  I — >1  is 
more  complicated,  with  Cj^q — >»for  e  -  0,  and  Cdq— *  0  for  e  >  0  due 
to  the  leading- edge  singularity.  For  e  =  0,  the  closed  model  has 


1 


vT^T 


Itl  _ 

- —  (\7  -  N/m> 

vf  -  l 


The  general  character  of  these  P-  and  Q-parameters  as  functions 
of  e  and  f  can  be  seen  in  Fig.  6  and  7  of  Ref.  1,  where  they  are  pre¬ 
sented  in  terms  of  the  flat  plate's  hydrodynamic  coefficients  for  the 
two  models,  using 


CDo/aZ  --  CDp,  CDci/a2  =  CDq,  CLJa 


CLp,  and  CLc1/q 


C 


DISCUSSION 


It  is  interesting  to  see  the  character  of  the  differences  between  the 
force  coefficients  for  the  two  models,  by  comparing  their  forms  for 
e  -  0.  The  open-model  relations  for  e  ~  0  given  above  lead  to 

CD  CL  TT 

Q2  Q  2 

while  the  corresponding  closed- model  relations  produce 
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la  contrast  to  these  results  for  free*  stream  speed  linearisation, 
the  well-known  cavity- wall- speed  linearisation  always  gives  linear 
variation  with  K  (not  K/a)  of  the  hydrodynamic  coefficients  for  K— *0. 
The  conversion  from  the  free- stream- speed  linearization  to  the  other 
is  quite  simple,  as  discussed  in  Appendix  A  of  Ref.  3.  It  is  easy  to  see 
that  the  Cp  conversion  given  there  for  e  =  1  applies  for  all  e  and  for 
both  the  closed  and  open  models.  Thus  for  a  given  profile,  cavity 
length,  and  particular  cavity  termination  within  the  open,  partly  closed, 
and  closed  model  family, 

(34)  Cj)cw  =  (1  +  KCw)  CDjg»  Clcw  =  (1+  Kcw)CL{8»  Cmcw  =  (1  +  Kcw)  CM{8 
where 

1  Kfs 

(35)  1  -  -  =  - 

n/TTkcw  2 

Thus,  for  example,  for  the  vented  flat  plate  with  e  =  0,  the  closed- 
model  result  is 

I  _  1  1 

-  (\TT  -  ~1).  1  -  -======  =  - - 

sim  v  i  +  kcw  n/fti 

Figure  17  shows  that  the  consequent  differences  in  the  curves  of  Cj^/a^ 
and  Cj_,/q  versus  K  are  not  important  for  the  small  angles  of  attack 
typical  of  vented  hydrofoil  applications. 

Applications  of  the  six  flat-plate  parameters  are  illustrated  next 
by  generation  of  the  CD  and  Cl  versus  K  curves  for  two  hydrofoil  pro¬ 
files  that  have  been  tested  in  fully  vented  flow  (Ref.  5).  As  discussed 
in  Ref.  1,  the  P-parameters  are  used  in  two  ways:  (1)  to  shift  from 
the  open-model,  smooth- entry  condition  to  the  cusp- closure  condition, 
and  (2)  to  obtain  the  open-model  results  for  arbitrary  attack  angle  of 
an  arbitrary  profile,  starting  with  its  cusp- closure  parameters.  The 
Q-parameters  are  used  to  obtain  the  closed- model  results  for  arbitrary 
attack  angle  of  an  arbitrary  profile,  starting  with  its  cusp-closure 
parameters. 


CD 


cw 


CLc 


w 


=  Ml  +  Kcw) 


0.0  0.1  0.2  0.3  0.4  0.  *>  0.6  0.7  0.8 

K 

FIG.  17.  Comparison  of  Lift  Coefficient  and  Cavity  Length  Versus  Cavity 
Number  for  Vented  Flat  Plate,  e  =  0  and  Closed  Model,  According  to 
Free-Stream-Speed  and  Cavity- Wall- Speed  Linearization. 
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EXAMPLES  OF  THEORETICAL  Cl  AND  Cd 
VERSUS  K  AND  TWO  COMPARISONS 
OF  THEORY  AND  EXPERIMENT 


HYDROFOIL  A 

The  profile  of  Hydrofoil  A  has  the  NACA  0010  thickness  distribu¬ 
tion  and  the  uniform- load  camber  line  for  design  Cl  of  0.  Z.  Thus  the 
s-parameters  for  Hydrofoil  A,  e  =  0.  3,  are  simply  obtained  from 
previous  results,  using,  in  obvious  notation, 

U<:«‘  (Uc“)o010  +  0-2(“'s)a=1.0 

and  corresponding  forms  for  a8,  tWs,  and  Cl8.  All  other  combinations 
of  the  NACA  four-digit  thickness  distribution  and  the  uniform-load 
camber  line  for  e  =  0.  3  can  also  be  considered  with  the  theoretical  re¬ 
sults  in  this  report  because  the  thickness  and  camber  distributions  were 
treated  separately.  If  Hydrofoil  A  had  been  treated  from  the  start,  only 
other  thickness  ratios  each  for  design  Cl  of  twice  the  thickness  ratio 
could  be  considered. 

As  discussed  in  Ref.  1,  the  cusp- closure  angle  of  attack,  a+,  is  a 
natural  and  convenient  reference  angle  for  each  cavity  length.  The 
values  of  a+,  uc+  and  CjJ  are  obtained  by  combining  the  P  solution  and 
the  particular  profile's  s- solution  for  each  1  to  make  tw  =  0.  Figure  18 


FIG,  18.  Cusp- Closure  Parameters  for  Hydrofoil  A, 
e  a  0.3. 
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gives  the  cusp- closure  parameters  for  Hydrofoil  A,  e  =  0.  3,  obtained 
with 

0  =  tw8  +  (a+  -  Qs)twp 


uc  =  uCg  - 


—  uCp.  c£  =  Cl8 - -  CLj 

*wp  *wp 


Having  K+  =  2u£  and  a+  for  each  i,  then  as  discussed  in  Ref.  1,  all 
a  values  for  a  >  a+  correspond  to  proper  cavity  solutions  in  the  re¬ 
stricted  sense  that  the  closure  singularity  has  proper  sign.  This 
merely  means  that  cases  of  definite  free- streamline  crossing  have  been 
eliminated,  and  of  course  it  is  likely  that  experimental  cavity  conditions 
will  not  extend  down  to  the  negative  K  region  typically  included  in  the 
proper  solution  domain,  as  seen  below.  A  more  detailed  discussion 
of  solution  domains  is  given  in  Appendix  B. 


The  proper  solution  domains  for  Hydrofoil  A,  e  =  0.  3  for  «*<_!  <_1. 1, 
are  given  in  Fig.  19a  and  b  for  the  open  and  closed  models,  respectively. 
The  vertical  bars  give  typical  experimental  K  values  for  what  was 
labeled  fully  vented  flow  in  the  water- tunnel  tests  (Ref.  5).  Since  the 
corresponding  cavities  extended  out  of  view  or  were  at  least  two  or 
chord  lengths  long,  there  is  a  first  indication  of  conflict  between  the 
theory  here  for  unbounded  flow  and  the  tunnel  experiments,  because 
the  theoretical  cavity  lengths  for  the  experimental  K  and  a  values  are 
typically  less  than  1.  5  chord  lengths.. 

Since  the  open- model  K,  ct  range  is  farther  from  the  experimental 
conditions,  only  the  closed  model  is  considered  in  the  following  detailed 
presentation  of  theoretical  results  and  comparison  with  experiment. 
However,  the  open- model  curves  essentially  coincide  with  the  closed- 
model  curves  but  terminate  at  very  low  K  values. 


As  discussed  in  Ref.  1,  simple  combinations  of  the  cusp-closure 
parameters  and  the  Q-parameters  lead  to  the  closed  model  Cl,  Cd> 
and  K  versus  f  and  a.  Thus  Fig.  20  and  21  present  the  curves  (solid 
lines)  of  Cd  and  Cl  versus  K  for  various  constant  a  values  for  Hydro¬ 
foil  A,  e  -  0.  3.  The  dashed  straight  lines  in  Fig,  21  are  the  contours 
of  constant  f. 


The  underlying  simplicity  of  the  structure  of  Fig.  21  is  noteworthy. 
For  each  t,  as  a  increases  from  a+,  one  moves  toward  higher  Cl  and 
K  along  a  straight  line  since 

Cl  -  c£  +  (a  -  a+)CLQ.  K  =  Kf  +  2(a  -  a+)uCQ 

Also  equally  spaced  a  values  appear  at  equal  intervals  along  each  con¬ 
stant  (  line,  the  interval  depending  on  f.  The  f  contour  curves  have 
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G.  20.  Theoretical  Cd  Versus  K  for  Hydrofoil  A,  e  =  0.3,  for  Various 
Values,  and  Comparison  With  Experiment.  Closed  model. 


FIG.  21.  Theoretical  CL  Versus  K  for  Hydrofoil  A,  e  =  0.3,  for  Various  a  Values 
and  Comparison  With  Experiment.  Closed  model. 
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been  omitted  from  Fig.  20  to  avoid  crowding,  but  it  can  be  seen  that 
these  would  have  parabolic  shape,  and  if  \l  Cd  were  plotted  versus  K, 
the  f  contours  would  be  straight  lines. 


However,  greater  insight  into  the  complicated  structure  of  Fig.  20 
can  be  obtained  more  easily  by  plotting  \l Cd  versus  a  as  in  Fig.  22. 

The  I  contours  are  again  straight  lines  and  the  intersections  with  the 
K  contours  are  spaced  along  each  i  contour  at  intervals  proportional 
to  the  K  increments.  Note  that  there  are  two  contours  for  K  =  0  in  the 
proper  solution  domain  [\l  Cd  >  0).  The  first  is  that  for  f  =°°,  with 
a  >  a+.  The  second,  for  finite  i,  is  similar  to  the  type  of  solution  for 
finite  i  and  K  =  0  which  was  discussed  in  Appendix  B  of  Ref.  2.  Since 
the  defining  equation  is 

K  =  0  =  2  [u+  +  (a  -  ci+)ucq] 

so  that  _ 

\TC5  -  (a  -  a+)  n/"Cdq  -  -  (uc/ucq)  VCDq 

the  intersection  of  the  two  contours  is  at  f  -  *>and 


=  n/  2tt 


lim  (-uc/uc0) 

f— »uo  ~ 


This  intersection  has  been  estimated  by  extrapolation  in  Fig.  22. 

The  various  types  of  experimental  data  points  in  Fig.  20  andT21 
and  the  choked- flow  curve  in  Fig.  20  are  discussed  later. 


HYDROFOIL  B 

Because  this  modified  NACA  65?-015  profile  is  uncambered,  the 
theoretical  treatment  for  Hydrofoil  B,  e  -  0.  3  is  similar  to  that  for 
the  NACA  0010  profile.  Thus  only  the  final  closed-model  curves  of 
Cd  and  Cl  versus  K  and  a  are  given  in  Fig.  23  and  24,  along  with  the 
corresponding  experimental  fully  vented  data  from  Ref.  5. 


DISCUSSION 

The  agreement  between  theory  and  experiment  seen  in  Fig.  20  to 
24  is  poorer  than  has  been  found  in  previous  cases  of  both  cavitating 
profiles  for  e  •-  0,  0.5,  and  1.0  and  of  vented  profiles  for  e  ~  0  and  1. 

It  is  believed  that  any  major  errors  in  the  theoretical  calculations 
would  have  been  detected  by  the  various  checks  described  earlier. 
While  the  accuracy  of  the  planimeter  integrations  decreases  for  small 
l  due  to  the  greater  complexity  of  Fj ,  ...  F4,  the  major  causes  of  the 
differences  between  the  theory  for  unbounded  flow  and  the  experiments 
arc  thought  to  be  real- fluid  and  tunnel  blockage  effects,  as  discussed 
below. 
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Consider  first  the  dashed  curves  in  Fig.  20  and  23.  These  give  the 
exact  theoretical  choked- flow  relation  between  Cd>  K,  and  tunnel  height- 
to- chord  ratio,  t,  from  Ref.  6: 

(36)  CD  =  N  1  +  K  -  l)2  r 

For  the  tests  of  Ref.  5,  r  =  14/4  =  3.  5.  Since  the  open  data  points  in 
the  figures  are  apparently  for  choked  or  nearly  choked  flow,  as  dis¬ 
cussed  later,  the  large  departure  of  those  experimental  Cq  values  from 
the  dashed  curves  shows  that  important  discrepancies  are  present  be¬ 
tween  the  real  and  the  ideal  flow  conditions.  Boundary- layer  separa¬ 
tion  on  the  unvented  lower  side  appears  to  be  a  major  cause  of  the  in¬ 
crease  of  Cd  as  a  decreases  from  zero.  This  explanation  is  suggested 
by  the  similar  increase  of  Cd  for  fully  wetted  flow,  seen  in  Fig.  12  of 
Ref.  5,  where,  for  example,  Cd  for  a  =  -6.4  degrees  is  higher  by  about 
0.0084  than  it  is  for  a  =  0.  Furthermore,  with  venting  it  seems  rea¬ 
sonable  to  expect  more  forward  boundary- layer  separation  on  the  lower 
side,  and  comparison  of  the  lower- surface  pressure  distributions  with 
and  without  venting  in  Fig.  30  of  Ref.  5  seems  to  silpport  this  idea. 
Finally,  if  the  experimental  Cd  values  are  reduced  by  0.0059  to  allow 
for  skin  friction,  as  suggested  in  Ref.  5,  page  106,  the  remaining 
separation  of  the  Cd  values  for  large  a  from  the  choked- flow  curve 
may  be  due  to  the  effect  of  blockage  on  the  experimental  K  determina¬ 
tion,  as  discussed  in  Ref.  5  (page  104). 

There  are  several  reasons  for  describing  the  fully  vented  condi¬ 
tions  in  the  tests  of  Ref.  5  as  being  apparently  choked  flow.  First 
there  is  the  characteristic  behavior  of  K  versus  gas  supply  (Fig.  28  of 
Ref.  5).  More  conclusive,  however,  is  the  notably  linear  variation  of 
K,  Cl  and  Cm  versus  a  (or  a.  Cl  and  Cm  versus  K)  as  seen  in  Fig. 56, 

7,  and  13,  respectively  of  Ref.  5,  That  this  behavior  is  precisely  that 
of  choked  flow  can  be  seen  as  follows.  As  long  as  cavity  length  is 
constant--as  it  is  of  course  in  choked  flow  when  I  =  oo--then  as  a  is 
varied  by  adding  the  flat  plate  solution  for  Aq,  K,  Cl  and  Cm  will  vary 
linearly  with  a. 

A  third  possible  cause  of  discrepancy  between  theory  and  experi¬ 
ment  is  as  follows.  A  key  difference  between  the  experiments  of  Ref.  5 
and  previous  vented-flow  tests  is  the  exhausting  of  gas  at  an  arbitrary 
chordwise  position  on  the  profile  without  any  sharp  edge  from  which  the 
free  streamline  springs  clear  of  the  profile.  Thus  the  theoretical  gas- 
water  interference  is  often  actually  a  layer  of  individual  bubbles  near 
the  exhaust  point,  as  seen  in  Fig.  21  of  Ref.  5.  In  other  words,  par¬ 
ticularly  for  negative  a  and  exhaust  ahead  of  the  maximum  thickness, 
some  differences  are  to  be  expected  as  the  result  of  the  free- streamline 
separation  being  delayed  till  some  point  aft  of  the  exhaust  point. 

The  last  items  of  experimental  data  to  be  discussed  are  the  data 
points  with  x  inside  and  the  "tracks"  joining  the  clear  and  the  x  data 
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points  in  Fig.  20  and  21.  In  Ref.  5,  the  x  data  identify  partially  vented 
flow  based  on  visual  observation,  and  the  "tracks"  indicate  unsteady 
forces.  A  definite  question  about  the  visual  classification  of  cavity 
length  is  implied  by  the  agreement  of  the  slopes  of  the  full-cavity  theory 
and  the  "track"  data  up  to  K  ~  0.  3,  where  corresponding  jogs  in  the  ex¬ 
perimental  Cl  curves  and  slope  changes  in  the  Cq  curves  seem  to 
occur.  This  suggests  that  the  lower  K  portions  of  the  "track"  data, 
including  some  of  the  x  data  points,  correspond  to  full- cavity  flow. 

In  discussing  this  seeming  conflict  of  the  visual  and  theoretically 
implied  indications  of  flow  configuration  with  the  authors  of  Ref.  5,  the 
following  hypothesis  was  developed.  Figure  25  shows  the  typical  se¬ 
quence  of  flow  configurations  observed  in  the  tests  of  Ref.  5  for  in¬ 
creasing  and  decreasing  gas  supply.  It  is  believed  that  configurations 
B  and  E  can  have  appreciably  different  gas  supply  rates,  E  being  the 
lower,  and  yet  that  they  are  essentially  the  same  with  respect  to  cavity 
pressure  and  hydrodynamic  coefficients.  However,  they  would  be 
classified  visually  as  partially  vented  and  full-cavity  flow,  respec¬ 
tively.  Thus  the  hypothesis  is  that  the  lower  K  portions  of  the  "track" 
data  of  Ref.  5  effectively  correspond  to  full-cavity  flow,  despite  the 
visual  classification,  and  that  the  qualitative  agreement  with  the  full- 
cavity  theory  is  to  be  expected.  Related  behavior  has  occurred  in  other 
cases  of  clear  and  frothy  cavities,  for  example  as  reported  in  Ref.  8, 
where  the  theoretical  and  the  measured  drags  continue  to  agree  as  the 
cavity  becomes  frothy  and  then  becomes  essentially  a  separated  flow 
region  filled  with  bubbly  water,  as  reported  in  Ref.  9. 


A 

,  .. . 


Z\  D" 

FIG.  25.  Sequence  of  Experimental  Cavity  Configurations 
for  Increasing  and  Decreasing  Gas  Supply. 
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SUMMARY 

The  example  calculations,  comparisons  of  theory  with  experiment, 
and  function  tables  in  this  report  are  the  culmination  of  a  study  of  the 
linearized  theory  of  vented  hydrofoils  in  two-dimensional  flow. 

The  tables  simplify  the  lift  and  drag  calculations  for  nearly  arbi¬ 
trary  hydrofoil  profiles  and  typical  exhaust  locations  for  infinite  or 
finite  cavity  length  of  full  cavity  flow.  Moment  coefficient,  because  of 
its  greater  calculation  complexity,  is  treated  only  for  infinite  cavity 
length. 

Two  experimental  hydrofoils  have  been  treated  with  the  theory  to 
illustrate  the  calculation  procedure  and  the  nature  of  the  theoretical 
curves  of  Cl  and  Cd  versus  K  and  a  for  e  =  0.  3.  The  resultant  agree¬ 
ment  of  theory  and  experiment  is  poorer  than  expected  from  previous 
comparisons.  However,  the  main  discrepancies  are  reasonably  at¬ 
tributed  to  hitherto  less  important  real-fluid  and  tunnel- choking  effects. 
In  particular,  boundary- layer  separation  on  the  lower,  unvented  side 
seems  to  have  been  appreciable  for  negative  attack  angles,  possibly 
aggravated  by  the  upper  side  venting. 

While  no  definite  statement  about  the  importance  of  the  choked- flow 
effects  upon  the  present  comparisons  can  be  made  until  choked- flow- 
theory  calculations  are  completed,  the  experimental  data  display 
qualitative  features  of  choked  flow.  The  indication  of  greater  effects 
of  choking  than  noticed  before  seems  reasonable  in  view  of  the  lower 
tunnel  height- to- chord  ratio  and  the  lower  theoretical  cavity  lengths  for 
unbounded  flow  for  the  experimental  choked  flow  cavity  numbers. 

An  unexpected  result  of  the  comparison  for  Hydrofoil  A  is  a  new 
interpretation  of  the  visual  classification  of  the  experimental  flows  as 
fully  vented  or  partially  vented,  following  from  the  qualitative  agree¬ 
ment  of  portions  of  the  data  for  so-called  partially  vented  flow  with  full 
cavity  theory. 
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Appendix  A 
CORRIGENDA 


Since  this  report  completes  the  theoretical  work  on  vented  hydro¬ 
foils  of  Ref.  1,  2,  3,  and  10,  errors  that  have  been  noted  in  two  of  those 
reports  are  corrected  here. 


Item  1:  In  Ref.  2,  page  15,  the  coefficient  of  d2  in  the  Cm  formula 
should  be  2  (1  +  6  cos^  0')  instead  of  2  cos^0'(5  +  2  cos^0'). 

Item  2:  In  Ref.  2,  page  18,  each  of  the  group  of  five  formulas  for 
(a)Cjj  =  o,  labeled  Case  Ha  through  IVb,  should  be  re¬ 
versed  in  sign. 

Item  3:  In  Ref.  10,  page  5,  the  formula  for  -Cm/s  is  incorrect 
and  should  read 


CM 
Q  2 


ir  / 1  +  \fe \4 


f 

1  +  4 

1  1  -  \e ' 

2 

1 

1  +  sfe  , 

The  corresponding  curve  on  the  same  page  should  be 
change  accordingly.  The  corrected  curve  rises  slightly 
above  the  straight  line  that  would  connect  the  end  points, 
which  remain  unchanged,  rather  than  slightly  below  as 
given.  Also,  the  corresponding  flat-plate  term  in  the 
Cm  formula  for  an  arbitrary  profile  on  page  6  should  be 
changed  to  read 


CM  " 


I  1  -  *Je 

1+41 - 

\  1  +  \fe 


2 


Item  4:  In  Ref.  10,  page  7,  the  factor  2/tt  in  the  ^  Cq  formula 

should  be  \'2/tt  and  the  coefficient  of  _  h  (  - -)d0 

o  {°  v  1 ta  ' 

in  the  Cm  formula  should  be  a/(l  +  a)i  instead  of 
a(l  +  2aZ  -  2a4)/ (1  +  a)3. 


Suggested  by  C.  J.  Henry  of  Davidson  Laboratory,  Stevens  In¬ 
stitute  of  Technology. 
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Appendix  B 
VENTED  WEDGE 


The  vented  wedge  is  almost  as  simple  as  the  flat  plate  to  treat  for 
finite  i,  because  again  point  singularities  on  the  unit  circle  produce  the 
solution. 

From  work  reported  in  Ref.  2,  the  smooth* entry  solution  for  the 
vented  edge  for  f  =  oo  is  seen  to  be 


6  /  4  -  e 

ws(4)  =  -  loge' - 


i0' 


18' 


giving 

and 


4  -  e 

Rl.  [w8(^1)]  =  uCs  =  0 

6 

-Im.  [wg(ei0)J  =  6-  a,  -  —  (it  -  0'),  0  <  0  < 


=  -qs  =  ~(-0'),  0'  <  0  <  IT 
tr 

Thus  for  finite  I  the  open* model,  smooth- entry  solution  wg(£)  with 

ws(4»)  *  0  is 


ws(4)  -  - 


loge 


4  -  e 


iO* 


—  +  loge  - 


4  *  ■  e 


■10’ 


\4  -  e-i9  /  \4oo-e1 

By  taking 

4^  -  e’10'  c  -  a  +  i(d  +  \fl~- ~aJ)  [(c  -  a)^  +  (d  +  ^1  -  a^)2]  *■  ei4- 


.16' 


4oo-  ei6'  c  -  a  +  i(d  -  y)  1  -  a2)  [(c  -  a)2  +  (d  -  yf\  -  a*)2]  *  e^' 

where 


tan  ^  = 

there  is  obtained 


d  +  Vl  -  a2 


c  -  a 


d  -  J 1  -  a2 

tan  4*'  - - 

c  -  a 
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and 


-Im.  [  w8(ei6)]  =  6  -  as  =  -  (ir  -  Q'  +  ij>'  -  4<),  0  <  0  <  9' 

XT 

6 

=  -  as  =  —  (-6'  +  tjj'  -  4>).  9'  <  0  <  it 

TT 


R1-  l  Ws^l>]  =  ucs  *  —  lo8 


2tt 


1  +  (tan  i}*)^  ] 
1  +  (tan  sj/)^ 


Proceeding  as  usual,  with  r,oo  -  (a  +  ib)  =  Ne*v,  it  is  found  that 

d£|  6tsfl  -  a2 


CLs  dW 

Vws  +  1 -  =  2ir - 

8  2  dr, 


dZ 


’ao 


N 


(sin  v  +•  i  cos  v) 


Though  the  mapping  for  <  >  1  is  singular  for  I — *1,  the  s-parameters 
are  not.  Ordinary  procedures  show  that  the  limit  values  for  J  =  1  are 


TTU, 


1 


—  =  -  loge 
6  2 


1  + 


1  +  \I~e  I  2  \fe 
1  -  \fe  \  1  -  \Te 


X 

1  +  \fe 

2  *fe 

2 

y  1  -  sfe 

\  1  -  \Te 

iras 


tan 


-1 


1  +  */e 
1  -  s /  e 

^s  I 

-  =  e4 


*  J— 


-  tan 


-1 


/ 1  -  e  I 

!1  +  's/e  Cr 


|1  +  *s/e 

1  -  \le 


1  ~  M  e 


=  e* 


> 2  \Je 
1  -  \Fe  , 


The  s-parameter  integrals  for  the  vented  wedge,  with  v  =  6  on  the 
upper  wetted  surface  and  v  =  0  on  the  lower,  are  easily  seen  to  be 


n-uCg  +  i-rras  =  6 


r 

Jo 


(Ft  +  i  Fz)  d0 


CL  bi 

3  2  N  JO 


f 


(F3  +  i  F4)  d0 


Figure  26  compares  the  analytic  and  planimeter  results  for  e  =  0.  5 
Thus  the  integral  forms  are  checked  and  the  over- all  accuracy  of  the 
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FIG.  26.  Analytically  Calculated  s- Parameters 
Versus  l/l  for  Vented  Wedge,  for  e  =  0.5,  and 
Comparison  With  Planimeter  Results. 


planimeter  calculation  procedure  is  demonstrated.  As  usual,  the 
curves  to  be  integrated  were  plotted  on  8-  by  10.5-inch  paper  with 
0  <  8  <  180°  along  18  centimeters  of  the  longer  side.  Thus  the  wedge 
check  on  accuracy  is  believed  pessimistic,  since  H  =  0  for  0  >  0'  as  60°. 
The  loss  of  accuracy  for  l  -  1.1  is  due  to  the  fact  that  nearly  equal 
areas  of  opposite  sign  happen  to  be  involved.  Because  fairing  of  the 
curves  was  a  much  greater  cause  of  error  than  the  planimeter  itself, 
the  additional  0  values  that  typically  had  to  be  considered  in  treating 
the  more  complex  H  or  H  functions  also  probably  increased  the  accu¬ 
racy.  However,  for  future  work,  a  basic  interval  of  5  degrees  would 
be  recommended. 

Because  the  wedge  solution  for  the  full  range  I  <  l  <  oo  has  been 
reached,  it  serves  as  a  useful  example  for  considering  the  coverage 
in  K  and  a  of  the  proper  open  and  closed  models,  according  to  the  cri¬ 
terion  from  Ref.  1  that  a  >  a+.  Thus  Fig.  27  gives  the  cusp-closure 
parameters  for  e  =  0.  5,  and  Fig.  28  shows  how  for  each  I  the  proper 
solutions  cover  the  line 
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K 

6 


K+ 

6 


a  -  a+ 

—  )  U*P.  Q 


The  open- model's  multivaluedness  for  the  typically  unimportant  area 
of  a  near  a"*"  means  that  a  given  K  and  a  can  be  obtained  with  more  than 
one  value  of  f,  for  which  the  hydrodynamic  coefficients  will  also  differ 
somewhat.  Actually  the  criterion  for  solution  admissibility  of  a  >  a+ 
is  not  a  clear-cut  condition  since  the  implied  free  streamlines  should 
be  calculated  and  any  cases  where  the  free  streamlines  cross  even  for 
a  slightly  greater  than  a+  should  be  excluded.  That  such  cases  occur 
is  suggested  by  considering  the  shape  of  the  free  streamlines  for  cusp- 
closure,  as  follows.  Since  the  profile- plus- cavity  obstacle  keeps  the 
same  pressure  distribution  if  the  free  stream  is  considered  reversed, 
it  would  appear  that  for  u^  positive,  as  sometimes  occurs,  the  free 
streamlines  should  initially  diverge  from  the  cusp  in  a  crossing 
manner.  However,  such  details  have  not  been  found  important  in  the 
practical  application  of  the  theory. 


In  a  more  practical  respect,  only  solutions  for  K  >  0  are  of  interest 
for  cavitating  flows  and  even  then  only  if  pressure  coefficient  on  the  foil 
nowhere  falls  below  -K.  For  vented  flows,  for  which  cases  of  negative 
K  can  be  conjectured,  it  remains  to  be  seen  whether  real  fluid  effects 
will  not  prevent  such  cases  from  occurring,  even  if  both  the  linearized 
theory  and  corresponding  nonlinearized  potential- flow  solutions  are 
spatially  reasonable,  as  in  the  case  of  some  cusp- closure  cavities. 
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Appendix  C 

PARTLY  CLOSED  MODEL  SOLUTIONS  FOR 
VENTED  FLAT  PLATE 

As  discussed  in  Re£.  1,  for  each  vented  hydrofoil  profile  there 
exists  a  family  of  partly  closed  cavity  model  solutions  for  given  cavity 
length,  with  wake  thickness  ty,  lying  between  the  zero  of  the  closed 
model  and  tWo  of  the  open  model.  The  particular  case  of  the  vented 
flat  plate  with  e  =  0  is  considered  here  in  order  to  compare  the  theo¬ 
retical  ranges  of  cavity  length  and  wake  thickness  with  certain  experi¬ 
mental  results  given  in  Ref.  11. 

The  partly  closed  model  solution  for  unit  angle  of  attack  for  each 
particular  cavity  length  is  constructed  by  superposing  the  open-  and 
closed- model  solutions  wp  and  wq  with 

w  =  T  Wp  and  (1  -  T)  wq 

tw  =  T  twp,  uc  =  T  uCp  +  (1  -  T)  uCq 

and  with  the  z- plane  closure- singularity  strength 

Cz  =  0  -  T)N/CDQ/2Tr 

Consideration  of  the  implied  cavity  and  wake  geometry  and  the  drag 
integral  shows  that  the  drag  coefficient  for  the  partly  closed  model  is 

CD  =  ^  2uc  +  Ztt{Cz)^ 

-  T  tWp  l  [T  uCp  +  (1  -  T)uCq]  +  (1  -  T)^  CDq 

The  lift  superposition  is  simply 

CL  *  T  CLp  +  (1  -  T)CLq 

and  since  Cj^/a^  =  C^/n  for  the  flat  plate  with  e  =  0,  it  must  also  be 
(as  was  found  for  e  =  1  in  Ref.  3  and  was  stated  for  all  e  in  Ref.  1)  that 

CD  =  T  CDp+  (1  -  T)CDq 

This  is  confirmed  by  application  of  the  formulas  given  earlier  for  the 
open-  and  closed- model  solutions  for  e  =  0. 
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It  is  not  clear  from  the  resultant  partly  closed  model  relations, 
which  give  K  and  tyy  as  functions  of  i  and  T,  and  which  imply  U#  as  a 
function  of  T  and  K,  that  maximum  wake  thickness  for  given  K  cor¬ 
responds  to  the  open- model  limit.  However,  manipulation  of  the 
relations  yields 

dtwl  w  I  T  df  1  \ 

—  =-n/T  1+ - 


U  dT 


it  / 1  +  sir  - 1 

=  -  <s/T(l  -  T)  - 

2  \  f  -  1 


I  l  +  sTT  v/<TT  \  1  “ 1 

T  +  (1  -  T) - - - - -  I 


From  this  it  can  be  seen  that  tw  monontonically  increases  with  T  for 
constant  K,  reaching  a  maximum  and  critical  point  for  T  =  1. 

To  illustrate  the  partly  closed  model  results.  Fig.  29  gives  the 
curves  of  force  coefficient,  cavity  length,  and  wake  thickness  versus 
cavity  number  for  T  =  0,  1/2,  and  1.  In  Fig.  30  some  portions  of  the 
same  cavity- length  and  wake- thickness  curves  are  compared  with  some 
experimental  results  of  T.  E.  Dawson,  who  tested  a  flat  plate  for  e  =  0 
(Ref.  11).  The  upper  portion  of  Fig.  30  shows  that  the  experimental 
cavity  lengths  correspond  to  T  ~  3/4.  If  the  shortness  of  the  cavities 
is  due  to  the  wake  produced  by  the  required  gas  supply,  rather  than  to 
the  effects  of  gravity,  the  free  surface,  or  the  bottom  of  the  tunnel, 
then  the  experimental  gas  supply  rates  should  be  compatible,  at  least 
roughly,  with  the  wake  thicknesses  corresponding  to  T  ~3/4.  The 
lower  half  of  Fig.  30  shows  the  experimentally  required  gas- supply  co¬ 
efficients,  Cq,  which  are  apparently  from  the  same  tests  as  for  the  f 
measurements.  The  definition  of  Cq  in  Ref.  11  makes  it  correspond 
to  tyj/a ,  since 

Q 

cQ  = - 

Uoo  be  sin  a 

where 

Q  =  volume  of  gas  supplied  to  test  hydrofoil  per  unit  time  (reference 
pressure  and  temperature  assumed  very  nearly  equal  to  free- 
stream  pressure  at  hydrofoil  depth  and  water  temperature) 

b  =  hydrofoil  span 

c  =  hydrofoil  chord 

In  the  lower  half  of  Fig.  30,  the  experimental  Cq  values  are  seen  to 
correspond  more  nearly  to  T  s  1/4  than  to  3/4,  but  in  view  of  the  crude¬ 
ness  of  the  wake  model,  the  order- of- magnitude  agreement  seems 
significant.  Also,  since  the  real  wake  near  the  cavity  is  a  mixture  of 
gas  and  water  moving  downstream  relative  to  the  foil  at  a  speed  less 
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than  Uoo.  it  is  to  be  expected  that  the  wake  thickness  implied  by  gas 
supply  would  be  smaller  than  that  implied  by  the  hydrodynamic  effect 
of  the  wake. 

It  will  be  interesting  to  see  whether  future  experiments  show  other 
cases  of  peculiarly  short  cavities  coupled  with  gas- supply  coefficients 
of  the  order  of  tWo/4  ss  Cd/4K,  in  accord  with  the  partly  closed  cavity 
model  interpretation  of  Dawson's  results. 
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1  Convair  Hydrodynamics  Laboratory,  San  Diego 
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1  Convair  Scientific  Research  Laboratory,  San  Diego 
1  Cornell  Aeronautical  Laboratory,  Inc.,  Buffalo  (Department  410) 

1  Cornell  University,  Graduate  School  of  Aeronautical  Engineering, 
Ithaca  (Library) 

2  Davidson  Laboratory,  Stevens  Institute  of  Technology,  Hoboken,  N.  J. 

A.  Suarez  (1) 

Dr.  J.  Breslin  (1) 

1  Douglas  Aircraft  Company,  Inc.,  Long  Beach  (Aerodynamics  Section) 

1  Eastern  Research  Group,  New  York  City 
1  EDO  Corporation,  College  Point,  N.  Y.  (Library) 

1  Electric  Boat  Division,  General  Dynamics  Corporation,  Groton,  Conn. 
(Library) 

1  Engineering  Societies  Library,  New  York  City 

1  General  Electric  Company,  Defense  Electronics  Division,  Pittsfield, 
Mass.  (Engineering  Librarian) 

1  General  Electric  Company,  Schenectady  (Librarian,  LMEE  Depart¬ 
ment) 

2  Gibbs  and  Cox,  Inc.,  New  York  City 

Dr.  S.  Hoerner  (1) 

Library  (1) 

1  Grumman  Aircraft  Engineering  Corporation,  Bethpage,  N.  Y.  (Library) 
1  Harvard  University,  Cambridge,  Mass.  (Department  of  Engineering 
Sciences) 

1  Hughes  Aircraft  Company,  Culver  City,  Calif.  (Library) 

5  Hydrodynamics  Laboratory,  CIT,  Pasadena 
Prof.  A.  J.  Acosta  (1) 

Prof.  A.  T.  Ellis  (1) 

Prof.  T.  Y.  Wu  (1) 

T.  Kiceniuk  (1) 

Library  (1) 

1  Hydronautics,  Inc.,  Rockville,  Md. 

1  Illinois  Institute  of  Technology,  Chicago  (Head,  Department  of 
Mechanical  Engineering) 

1  Institute  of  Aerospace  Sciences,  Inc.,  New  York  City  (Librarian) 

1  Johns  Hopkins  University,  Baltimore  (Head,  Department  of  Mechani¬ 
cal  Engineering) 

1  Lehigh  University,  Bethlehem,  Pa.  (Civil  Engineering  Department) 

1  Lockheed  Aircraft  Corporation,  Burbank,  Calif.  (Library) 

1  Lockheed  Aircraft  Corporation,  Missiles  and  Space  Division,  Palo 
Alto,  Calif.  (R.  W.  Kermeen) 

2  Massachusetts  Institute  of  Technology,  Cambridge 

Department  of  Civil  Engineering  (1) 

Department  of  Naval  Architecture  and  Marine  Engineering  (1) 

1  McDonnell  Aircraft  Corporation,  St.  Louis  (Library) 

1  Mississippi  State  University,  State  College  (Aerophysics  Department) 

1  New  York  State  Maritime  College,  Fort  Schuyler  (Librarian) 

1  New  York  University,  Institute  of  Mathematical  Science,  New  York 
City 

1  North  American  Aviation,  Inc.,  Los  Angeles  (Technical  Library, 
Department  56) 
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1  Oceanic*  Incorporated,  New  York  City 

1  Operations  Research,  Inc.,  Los  Angeles 

2  Ordnance  Research  Laboratory,  Pennsylvania  State  University, 
University  Park  (Garfield  Thomas  Water  Tunnel) 

1  Pacific  Aeronautical  Library  of  the  IAS,  Los  Angeles 
1  Reed  Research,  Inc. 

1  Rensselaer  Polytechnic  Institute,  Troy,  N.  Y.  (Department  of  Mathe¬ 
matics) 

1  Republic  Aviation  Corporation,  Farmingdale,  N.  Y. 

1  Rocketdyne,  Canoga  Park,  Calif.  (Library,  Department  596-6) 

1  Rose  Polytechnic  Institute,  Terre  Haute,  Ind.  (Library) 

1  Society  of  Naval  Architects  and  Marine  Engineers,  New  York  City 

2  Southwest  Research  Institute,  San  Antonio 

Director,  Department  of  Applied  Mechanics  (1) 

Editor,  Applied  Mechanics  Review  (1) 

1  Stanford  University,  Stanford,  Calif.  (Department  of  Civil  Engineer¬ 
ing,  Prof.  B.  Perry) 

1  Technical  Research  Group,  Syosset,  N.  Y. 

1  The  Boeing  Company,  Seattle  (Library,  Organization  No.  2-5190) 

2  The  Martin  Company,  Baltimore  (Science  Technical  Library) 

1  The  Rand  Corporation,  Santa  Monica,  Calif.  (Technical  Library) 

4  The  University  of  Michigan,  Ann  Arbor 
Department  of  Civil  Engineering  (1) 

Department  of  Engineering  Mechanics  (1) 

Department  of  Naval  Architecture  and  Marine  Engineering  (2) 

1  The  University  of  Southern  California,  Los  Angeles  (Engineering 
Center) 

1  Thompson  Ramo  Wooldridge,  Inc.,  Cleveland  (Chief  Engineering 
Science  Group) 

1  United  Technology  Corporation,  Sunnyvale,  Calif.  (Technical  Library) 

2  University  of  California,  College  of  Engineering,  Berkeley 

Prof.  J.  V.  Wehausen  (1) 

Library  (1) 

1  University  of  Illinois,  Urbana  (College  of  Engineering) 

2  University  of  Iowa,  Iowa  Institute  of  Hydraulic  Research,  Iowa  City 

Prof.  H.  Rouse  (1) 

Prof.  L.  Landweber  (1) 

1  University  of  Kansas,  School  of  Architecture  and  Engineering, 
Lawrence 

1  University  of  Maryland,  Institute  of  Fluid  Dynamics  and  Applied 
Mathematics,  College  Park 

1  University  of  Minnesota,  St.  Anthony  Falls  Hydraulic  Laboratory, 
Minneapolis 

1  University  of  Nebraska,  Lincoln  (Department  of  Engineering 
Mechanics) 

1  University  of  Notre  Dame  (Department  of  Engineering  Mechanics) 

1  University  of  Wisconsin,  Mathematics  Research  Center,  Madison 
(Prof.  L.  M.  Milne- Thomson) 

1  Vitro  Corporation  of  America,  Silver  Spring  Laboratory  (Librarian) 
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